Autoresonant control of the many-electron dynamics in nonparabolic quantum wells 
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The optical response of nonparabolic quantum wells is dominated by a strong peak at the plas- 
mon frequency. When the electrons reach the anharmonic regions, resonant absorption becomes 
inefficient. This limitation is overcome by using a chirped laser pulse in the autoresonant regime. 
By direct simulations using the Wigner phase-space approach, we prove that, with a sequence of 
just a few pulses, electrons can be efficiently detrapped from a nonparabolic well. For an array of 
multiple quantum wells, we can create and control an electronic current by suitably applying an 
autoresonant laser pulse and a slowly varying dc electric field. 



Small semiconductor devices, such as quantum dots 
and quantum wells, have attracted considerable attention 
in recent years, particularly for possible applications in 
the emerging field of quantum computing jh]. For quan- 
tum devices working with many electrons [2] , it is crucial 
to understand the properties of the self-consistent elec- 
tron dynamics and its response to external electric fields. 
Particular attention has been devoted to intersubband 
transitions in semiconductor quantum wells, which take 
place on the meV energy scale and involve excitation fre- 
quencies of the order of the teraherz [3[ . Several theoreti- 
cal and computational studies have investigated the elec- 
tron response, mainly using Hartree-Fock semiconductor 
Bloch equations [3] or density functional theory (DFT) 
Q. For perfectly parabolic confinement, the electron re- 
sponse is dominated by the Kohn mode @, 0] , consisting 
of rigid oscillations of the electron gas at the effective 
plasmon frequency. For nonparabolic confinement, the 
Kohn mode still dominates the initial response. How- 
ever, when the electrons reach the anharmonic regions, 
the resonance condition is lost and absorption becomes 
inefficient. 

In this Letter, we show that this limitation can 
be overcome by resorting to autoresonant excitation 
[1, Basically, autoresonant excitation occurs 

when a classical nonlinear oscillator is excited by an os- 
cillating force with slowly varying frequency: F{t) = 
ecos [uJo{t — to) -t- ia(i — to)^] : where e is the excitation 
amplitude and loq the frequency, which matches the lin- 
earized oscillator frequency; a is the rate of variation of 
the excitation frequency. For \a\ <C Wg and e above a 
certain threshold, the instantaneous oscillator frequency 
becomes "locked" to the instantaneous excitation fre- 
quency, so that the resonance condition is always sat- 
isfied. In that case, the amplitude of the oscillations 
grows indefinitely and without saturation, until of course 
some other effect kicks in. For the single-particle case the 
threshold behaves as eth ~ la]^/"^ [8|. 

In order to study the self-consistent electron dynamics, 
we make use of the Wigner representation of quantum 
mechanics. A mixed quantum state is represented by a 
function of the phase space variables plus time: f{x,v,t) 
(we deal with one-dimensional problems), which evolves 



according to the Wigner equation 
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where m^tc is the effective electron mass and V{x, t) is 
the total potential acting on the electrons. The latter 
is composed of three terms: (i) the confining poten- 
tial V"conf(a;); (ii) the autoresonant oscillating potential 
Kuto = xF{t); and (iii) the Hartree potential VH{x,t), 
which obeys Poisson's equation = (e^/e) J_^fdv, 
where e is the absolute electron charge and e is the effec- 
tive dielectric constant. We consider wide quantum wells 
(w lOOnm) at moderate electron temperatures (« 2Tf), 
for which the exchange and correlation corrections can 
be neglected [11, jl2] . 

The right-hand side of Eq. ^ models disorder or 
phonon scattering in the form of a friction-diffusion term 
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where 7 is the relaxation rate (inverse of the relaxation 
time Ti), and D^, are diffusion coefficients in ve- 
locity and real space respectively, which are related to 
the decoherence time T2. In order for Eq. ([1]) to pre- 
serve the positivity of the density matrix associated to 
the Wigner distribution function, the scattering term 
must be in Lindblad form [3]. This is automatically 
achieved [Ti'l if the above coefficients respect the inequal- 
ity D^D.^ > 72ft2/4m|. 

We focus on confining potentials that can be ap- 
proximated by a parabola at the bottom of the well: 
Konf (a;) — ^uiqItl^x'^ -)-..., where the frequency ujq can 
be related to a fictitious homogeneous positive charge of 
density uq via the relation u>q = e^no/msice. We then 
normalize time to coq^', space to the harmonic oscillator 
length Lho = y^Jrn^too', velocity to y^hujQ/m^; energy 
to hioo] and the electron density to ng. 

As initial condition, we take a Maxwell-Boltzmann 
velocity distribution with temperature > Tp and a 
Gaussian density profile with peak density rig. This is 
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not an exact stationary state, but it evolves very little 
if no perturbation is applied. (The precise form of the 
initial state is irrelevant for our purposes, provided it is 
localized at the bottom of the well.) We define the filling 
fraction as rj = Ue/no < 1: the limit case 77 = cor- 
responds to very dilute densities, for which the Hartrec 
potential is negligible. 

We use typical parameters for semiconductor quan- 
tum wells 0: effective electron mass and dielectric con- 
stant = 0.067me and e = ISsq; volume density 
no = 5 X lOi^cm-3, luq = 1.35 x lO^^g-i^ ^ 8.9meV, 
Lho = ll.Snm. For 77 = 1, this yields a maximum sur- 
face density for the electrons Us — 1.35 x lO^^cm"^ and 
a maximum Fermi temperature Tp — 85. 7K. The elec- 
tron temperature is taken to be Te = 2hujo — 200K. The 
relaxation time is Ti « 70ps, which corresponds to a re- 
laxation rate 'j/ujq — 0.001. The diffusion coefficient in 
velocity space is Dy — ^ ^Jli^T^Jm^ . 

First, we consider a single quantum well with a Gaus- 
sian potential: 14onf(a^) = — exp(— x^/2a-^), with Vo — 
a'^m^ifUj^, a = 4Lho — 45nm, and an overall width of 
the quantum well of L = 24Lho — 270nm. This is in line 
with recent experiments on wide parabolic quantum wells 
(WPQWs), which can reach a width of several hundred 
nanometers [T2| . 

The electron gas is excited with a chirped pulse, with 
a = -O.OOlii^g, to = bQ{)ujQ^, and amplitude e = 0.2 > 
Eth: in normalized units. The latter corresponds to an 
electric field of the order of O.lmV/nm, which can be 
easily achieved experimentally. Chirped pulses were also 
suggested to excite transitions in two-level quantum sys- 
tems : a situation quite different from the wide quan- 
tum wells considered here, where many levels are present 
and progressively excited by means of the autoresonant 
technique. 
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FIG. 1: Single Gaussian well. Total number of electrons in 
the well, normalized to its value at f = 0, for different values 
of the filling fraction: 77 = (solid line), 77 = 0.1 (dotted), 
r] = 0.5 (dashed), rj — 1 (dot-dash). Inset: Time evolution of 
the dipole for 77 — 0.5. 

Figure 1 shows the dipole of the electron distribution 
(average of x) and the total number of electrons in the 



well (integral of f{x,v,t) in the phase space) normal- 
ized to its initial value, for different values of the fill- 
ing fraction 77. The dipole starts increasing rapidly at 
t = to, i.e. when the time-dependent frequency of the 
laser field crosses the linear resonance uiq. At cuqI ~ 900, 
the electrons that have been accelerated by the laser field 
are finally ejected from the well, leaving behind a rem- 
nant population that is still sitting at the bottom of the 
well. The dipole then decreases rapidly and remains close 
to zero thereafter. The number of electrons in the well 
is reduced of one order of magnitude at the end of the 
run. The slight initial decrease observed for large elec- 
tron densities {rj = 1) represents a small evaporation due 
to Coulomb repulsion. 
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FIG. 2: Time evolution of the total number of electrons for a 
series of four laser pulses; rj = 0.5. 

By applying a sequence of four identical pulses it is 
possible to reduce the number of electrons by almost three 
orders of magnitude, as shown in Fig. 2. We stress that, 
for a non-chirped pulse (a = 0) or for a chirped pulse 
below the autoresonant threshold (e < eth), virtually no 
electrons leave the well. In addition, the effect is still 
observed when the laser frequency is mismatched with 
respect to the harmonic oscillator frequency loq. Even 
for a mismatch of ±10%, the same number of electrons 
are ejected from the well. 

As a second application, we consider a periodic array 
of quantum wells, with a cosinusoidal confining poten- 
tial Vconf — — Vb cos(27ra;/A), where A is the width of 
each well, and Vq ~ mj^uj'^X'^ /2tt. Such periodic superlat- 
tices can be practically realized as multilayer semiconduc- 
tor heterostructures and there have been recent at- 
tempts at simulating these structures using Bose-Einstein 
condensates trapped in optical lattices [ij]. Here, we 
neglect the Hartree potential and consider noninteract- 
ing electrons (this amounts to assuming 77 <C 1). We 
take 7 = 0.002wo, a = -0.002wg, tuoto = 300, and 
e = 0.2fia;o/iho- Initially, each quantum well is occu- 
pied by a single electron, represented by a minimum un- 
certainty packet. We use a periodic computational box 
with spatial period equal to 3A. 

The idea is to create an electron current by applying 
a suitable laser pulse. From the previous study, we have 
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FIG. 3: Evolution of the average velocity for an array of quan- 
tum wells with dc bias. Inset: Velocity distribution of the 
electron population at ujo t = 0, 500 , 1000, and 2500. Velocity 
is measured in units of i/fttJo/wi*. 

learnt that the electrons can be efficiently extracted from 
the well. By applying a small constant (dc) electric field, 
we can force the electrons to always leave the well from 
one side and thus create an electric current. For this, we 
use a small dc field {Eq — O.lftwo/eLho — 0.08mV/nm) 
to polarize the array of quantum wells [ll|. In Fig. 3, we 
plot the evolution of the average velocity of the electron 
distribution {v) — J fvdvdx/ J fdvdx. 




FIG. 4: Evolution of the average velocity for an array of quan- 
tum wells with dc bias, in the case of two laser pulses. The 
dashed line represents the evolution of the dc field, in arbi- 
trary units. 



At t — to, Si, positive current starts building up and, 
by the end of the run, the average velocity reaches a 
value (v) ~ 25 in our normalized units. This value cor- 
responds to the asymptotic velocity reached by an elec- 
tron subjected to a constant electric field Eq = 0.1 in 
the presence of a friction coefficient 2j = 0.004, i.e. 
i^max = Eo/2^ — 25 (normalized units). This shows that 
all electrons have been detrapped from the potential well 
and accelerated by the dc field. The electron velocity 
distribution (Fig. 3, inset) conffi'ms this scenario. 

In order to make the whole process reversible, one 
would need to re-trap the electrons inside the potential 
well. This can be achieved in the following way (Fig. 
4): (i) after the current has been generated, the oscil- 
lating pulse is switched off suddenly at uj^t = 2000; (ii) 
at the same time, the dc electric field Eq is switched off 
adiabatically between loqI = 2000 and luqI — 4000 (expo- 
nential decrease with time constant r = 400^^^): during 
this phase, all electrons are re-trapped and the current 
goes back to zero; (iii) at loqI — 4000, the dc electric 
field is switched on again and another autoresonant laser 
pulse (identical to the first one) is used to excite the cur- 
rent once again. This procedure can be repeated several 
times, so that the electric current can be switched on and 
off, and can even change sign if one changes the sign of 
the dc field. 

In summary, we presented numerical evidence that the 
electron gas in wide nonparabolic quantum wells can be 
efficiently excited with an autoresonant laser pulse. Non- 
linear effects were triggered using a relatively weak pulse, 
even when the laser frequency is poorly tuned. These 
techniques could be used to achieve better control of the 
electron dynamics in quantum solid-state devices. 
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